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ABSTRACT 

We report the first submillimeter interferometric observations of an ultraluminous infrared 
galaxy. We observed Arp 220 in the CO J=3-2 line and 342 GHz continuum with the single 
baseline CSO-JCMT interferometer consisting of the Caltech Submillimeter Observatory (CSO) 
and the James Clerk Maxwell Telescope (JCMT). Models were fit to the measured visibilities to 
constrain the structure of the source. The morphologies of the CO J=3-2 line and 342 GHz 
continuum emission are similar to those seen in published maps at 230 and 110 GHz. We 
clearly detect a binary source separated by ~ 1" in the east- west direction in the 342 GHz 
continuum. The CO J=3-2 visibility amplitudes, however, indicate a more complicated structure, 
with evidence for a compact binary at some velocities and rather more extended structure at 
others. Less than 30% of the total CO J=3-2 emission is detected by the interferometer, which 
implies the presence of significant quantities of extended gas. We also obtained single-dish CO 
J=2-l, CO J=3-2 and HCN J=4-3 spectra. The HCN J=4-3 spectrum, unlike the CO spectra, 
is dominated by a single redshifted peak. The HCN J=4-3/CO J=3-2, HCN J=4-3/HCN J=l-0 
and CO J=3-2/2-l line ratios are larger in the redshifted (eastern) source, which suggests that 
the two sources may have different physical conditions. This result might be explained by the 
presence of an intense starburst that has begun to deplete or disperse the densest gas in the 
western source, while the eastern source harbors undispersed high density gas. 

Subject headings: galaxies: ISM - galaxies: indiviSual: (Arp 220) - galaxies: starburst - techniques: 
sub- mm interferometry 



1. Introduction 

Ultraluminous infrared galaxies contain ex- 
traordinary nuclear starbursts and, at least in 
some cases, an active galactic nucleus, all hid- 
den within a dense shroud of gas and dust. Un- 
derstanding the mechanisms which produce and 
power these luminous galaxies has taken on added 
urgency with the discovery of their young coun- 
terparts at cosmological distances (Ivison et al. 
2000). The nearest and prototype ultraluminous 
infrared galaxy, Arp 220, is located at a distance 
of 73 Mpc (H a = 75 km s _1 ) and is one of the 
best studied of this class of galaxies. The pres- 
ence of tidal tails observed in the optical (Arp 
1966; Joseph & Wright 1995) as well as two com- 
pact emission peaks at near-infrared (Scoville et 
al. 1998), millimeter (Scoville, Yun, & Bryant 
1997, hereafter SYB; Downes & Solomon 1998, 
hereafter D&S; Sakamoto et al. 1999, hereafter 
S99), and radio wavelengths (Bccklin & Wynn- 
Williams 1987; Norris 1988; Sopp & Alexander 
1991) suggests that Arp 220 is a recent merger. 
Given the observed correlation between high in- 
frared luminosity and disturbed optical morpholo- 
gies indicative of galaxy interactions (Mirabel & 
Sanders 1989), this merger is likely responsible for 
the extremely high infrared luminosity of Arp 220 
[1.4 x 10 12 L , Soifer et al. (1987)]. There has been 
much debate as to whether the high infrared lu- 
minosity is due to a starburst or an active nucleus 
or a combination of the two, both for Arp 220 in 
particular and for ultraluminous infrared galaxies 
in general (Genzel et al. 1998; SYB; Lutz et al. 
1996). Recent radio studies of Arp 220 provide 
support for the starburst hypothesis: the 18 cm 
flux seen in VLBI observations is thought to be 
emitted by luminous radio supernovae (Smith et 
al. 1998) and the effects of these supernovae winds 
are seen in X-rays (Hcckman et al. 1996). Near- 
infrared (Emerson et al. 1984; Rieke et al. 1985; 
Sturm et al. 1996; Lutz et al. 1996) and earlier 
radio observations (Condon et al. 1991; Sopp & 
Alexander 1991; Baan & Haschick 1995) further 
support the starburst scenario. 

Besides its extremely high infrared luminosity, 
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Arp 220 also contains large amounts of dust and 
molecular gas (~ 9 x 10 9 M Q SYB); large gas 
masses of 4 — 40 x 10 9 M Q are typical for ultra- 
luminous infrared galaxies (Sanders et al. 1991). 
The extinction at optical wavelengths is estimated 
to be at least Ay ~ 50 and possibly as high as 
A v ~ 1000 (Sturm et al. 1996; D&S); even at near- 
infrared wavelengths (2.2/um), dust lanes obscure 
the possible nuclei (Scoville et al. 1998) These 
high extinctions mean that radio observations are 
needed to probe the deep interior regions. D&S 
have mapped Arp 220 in CO J=2-l and 1.3 mm 
continuum; besides detecting two emission peaks, 
they also see an extended disk (with full-width 
half-maximum extent of 2" x 1.6"). They inter- 
pret the emission peaks as the nuclei of the merg- 
ing galaxies embedded in a more extended disk of 
molecular gas. S99 refine the model further us- 
ing their CO J=2-l and continuum observations; 
in this model, the nuclei are each embedded in 
their own gas disk which is counter-rotating in the 
larger common disk. An alternative interpretation 
of the emission peaks as being due to crowding in 
the orbits of molecular gas and stars is given by 
Eckart & Downes (2001). 

This paper presents the very first interferomet- 
ric submillimeter observations of Arp 220 in the 
CO J=3-2 line and 0.88 mm continuum. In fact, it 
is the first paper to present data from submillime- 
ter interferometry of any extragalactic source. We 
also present single dish observations of Arp 220 of 
HCN J=4-3, a high density tracer. Furthermore, 
to determine the total CO flux of Arp 220, which 
could be partly resolved out by the interferome- 
ter, complementary single dish data were taken in 
CO J=3-2 as well as in the CO J=2-l line. These 
submillimeter observations allow us to penetrate 
deep into the interior of the nuclei, while at the 
same time tracing hotter gas, which may be more 
closely associated with the source of the ultralu- 
minous infrared luminosity. We describe the ob- 
servations in § 2 and the data analysis in § 3. The 
interferometric data were obtained with the single 
baseline CSO-JCMT interferometer, the only cur- 
rently available submillimeter interferometer with 
sufficient bandwidth to observe the broad emission 
lines of Arp 220. With data from a single baseline, 
mapping is not possible and so we analyze the data 
by making fits in the visibility plane. The data are 
compared to published data at lower frequencies 
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and to single dish data in § 4 and the conclusions 
are presented in § 5. 

2. Observations and Data Reduction 
2.1. Interfero metric Data 

The intcrfcrometric measurements were obtained 
with the CSO-JCMT interferometer on Mauna 
Kea, Hawaii, consisting of the Caltech Submil- 
limctcr Observatory (CSO) and the James Clerk 
Maxwell Telescope (JCMT) 2 . Usually these two 
telescopes are operated independently for single 
dish observations, but occasionally they are linked 
together as a submillimctcr interferometer with a 
single baseline of 164 m and a minimum fringe 
spacing of 1.1" at 350 GHz (Lay et al. 1994a; Lay 
1994b; Lay, Carlstrom, & Hills 1995, 1997). The 
CSO-JCMT interferometer is currently the only 
submillimctcr interferometer with a bandwidth 
broad enough to accommodate the <~ 800 km s _1 
wide line of Arp 220. 

The 12 CO J=3-2 transition and associated con- 
tinuum emission at 342.5 GHz of Arp 220 were 
observed on 1997 May 10 in good weather with 1.4 
mm of precipitable water vapor. The coordinates 
of Arp 220 used were a(2000)=15 /l 34 m 57/19, 
(5(2000)^23° 30' ll'/3 (originally a(1950)=15' 1 
32 m 46 sglj <5(1950)=23° 40' 07'.'9). To moni- 
tor the gain of the system, these observations 
were interleaved with those of the quasar 3C 345. 
Continuum and line observations of Arp 220 and 
3C 345 were alternated in time. We also observed 
the quasar 3C 273 in both single dish and inter- 
ferometric mode to determine the absolute flux 
calibration; the interferometric observations were 
additionally used to determine the shape of the 
passband. Figure 1 shows the (u, v) track of the 
CO J=3-2 data with a fringe spacing ranging from 
l'/2 to 5'/8. 

Both the Arp 220 and the 3C 345 data were 
divided by a passband created from observations 
of 3C 273. Changes in the relative gain of the 
interferometer as a function of time were then cal- 
ibrated by applying the gain curve derived from 
observations of the point source 3C 345 to the 

2 The JCMT is operated by the Joint Astronomy Centre in 
Hilo, Hawaii on behalf of the parent organizations Particle 
Physics and Astronomy Research Council in the United 
Kingdom, the National Research Council of Canada and 
The Netherlands Organization for Scientific Research. 
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Fig. 1. — Track of CO J=3-2 interferometric ob- 
servations in (u, v) coordinates. 

Arp 220 data. Each set of ten 10-second inte- 
grations was vector-averaged to produce a series 
of 100-second integrations with higher signal-to- 
noise ratios. 

The spectra displayed in Figure 2 were further 
averaged to 1000 seconds. (The sixth spectrum 
only contains 500 seconds of data.) Since the aver- 
age phase might drift over time scales longer than 
100 seconds due to instrumental or atmospheric 
phase shifts, the average phase over the whole fre- 
quency range of each 100-second sample was mea- 
sured and then subtracted from the phases in each 
of the individual velocity channels. The resultant 
line data were then averaged by deriving the arith- 
metic means of the sine and cosine components in 
each individual channel. 

By comparing single dish measurements of 
Mars with those of 3C 273, we obtained a flux 
of 12 ± 0.5 Jy for 3C 273, which agrees very well 
with the value of 12.2 Jy measured at 850 /im with 
SCUBA around the same time (Robson, Stevens, 
& Jenness 2001). We used the 3C 345 measure- 
ments immediately before and after the 3C 273 
data and obtained a 3C 345 flux of 0.90 Jy aver- 
aged over the emission at 342.5 and 339.5 GHz. 
After averaging and gain calibrating the Arp 220 
and 3C 345 data in exactly the same way, the 
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Fig. 2. — Intcrfcrometric CO J=3-2 spectra. The 
histogram shows the flux on a scale of to 2.4 Jy. 
The dashes are the phase, where the upper half of 
each plot covers the phase from — tt to tt. The ve- 
locities are relative to 339.5 GHz corresponding to 
a velocity of 5450 km s _1 in the radio definition. 
The first spectrum on the top left was taken at 
-2.3 h (HA) at (u, v) coordinates (-140 kA,90 kA), 
the top right at HA=-0.8 h (-170 kA,60 kA), 
middle left at HA=2 h (-160 kA,10 kA), middle 
right at HA=3.2 h (-125 kA,-10 kA), bottom left 
at HA=4.5 h (-70 kA,-20 kA), bottom right at 
HA=5.5 h (-25 kA,-25 kA). The continuum - mea- 
sured at slightly different (u, v) positions - was 
extrapolated to the above (u, v) positions (using 
our best fitting binary model) and the expected 
continuum emissions of 0.16, 0.31 0.28, 0.15, 0.2 
and 0.35 Jy have been subtracted from the respec- 
tive line data. 



scalar average of the 3C 345 data was deter- 
mined and from it the conversion factor from K 
to Jy. This conversion may vary between the 
continuum and CO data and with velocity bin be- 
cause the gain calibration can vary with frequency 
range. We derived conversion factors ranging from 
48 Jy K _1 to 59 Jy K _1 . The overall uncertainty 
in the calibration due to gain variations and un- 
certainties in the planet flux is estimated to be 
20% (Lay 1994b). 

Since the continuum data were measured at dif- 
ferent (u, v) positions than the line data, we ex- 
trapolated the continuum data using the binary 
model discussed in §3.3 and subtracted 0.16, 0.31, 
0.28, 0.15, 0.2 and 0.45 Jy from the six CO mea- 
surements (with increasing hour angle), respec- 
tively. All velocities quoted use the radio defini- 
tion (v ra d = c(v rest - v obs )jv rest ) with respect to 
the local standard of rest (lsr) . 

2.2. Single-Dish Data 

On 1997 July 17 and 18 we used the JCMT to 
observe emission lines from 12 CO J=3-2, 12 CO 
J=2-l, and HCN J— 4-3 under good weather con- 
ditions. The observations were made with the fa- 
cility 230 GHz (A2) and 350 GHz (B3) receivers, 
which had system temperatures of 270-360 K and 
400-600 K, respectively, in the center of the spec- 
tral window. The spectra were obtained with a 
chopping secondary mirror using a switch cycle of 
1 Hz and a beam throw of 60" in azimuth. We 
configured the DAS to give a spectral resolution 
of 1.25 MHz with a bandwidth of 1.86 GHz for B3 
and 0.95 GHz for A2. The relative errors in point- 
ing were small, ~ 2" rms for A2 and ~ 3" rms for 
B3, with significant systematic errors in the point- 
ing apparent when observing at elevations greater 
than 80°. Spectra obtained at such high elevations 
were not used in the subsequent analysis. 

Our data were flux calibrated through observa- 
tions of Mars and Uranus. We found significant 
differences between our derived aperture efficien- 
cies and the typical aperture efficiencies for the 
JCMT. We derive aperture efficiencies of 0.48 at 
230 GHz, 0.42 at 265 GHz, and 0.43 at 345 GHz, 
compared to the standard telescope values of 0.69 
at 230 GHz and 0.58 at 345 GHz. We adopt 20% 
as the error in the absolute calibration of our single 
dish observations. We present all our single dish 
data on the T^ temperature scale, which is most 
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appropriate for the compact nuclear emission seen 
in Arp 220 (SYB; S99; D&S). 
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Fig. 3.— CO 3=2-1 and CO J=3-2 spectra with 
a resolution of 22". For CO J=2-l the integrated 
intensity is 53 K km s _1 or 1730 Jy km s _1 . For 
CO J=3-2 the integrated intensity over the same 
22" beam is 45 K km s" 1 or 3700 Jy km s _1 . 

We used the package SPECX (Padman 1993) 
to reduce the spectra and used first-order base- 
lines in all cases. The CO J =2-1 spectrum is 
shown in Figure 3. We derive an integrated in- 
tensity of 53 K km s -1 , which for a conversion 
factor of 32.6 Jy K _1 (T^) corresponds to a flux of 
1730 Jy km s _1 . This value is significantly larger 
than the flux determined by Radford, Solomon, & 
Downes (1991) in a 12" beam (1040 Jy km s" 1 ). 
The interferometric maps of SYB and D&S show 
very compact CO J=2-l emission with a total ex- 
tent of roughly 4 x 4". The larger flux detected 
in the 21" beam of the JCMT suggests that there 
may also be extended CO J=2-l emission with 
roughly half the total intensity seen in the nuclear 
region. However, we also obtained CO J=2-l spec- 
tra at positions offset 11" (half a beam width) from 
the central position; the flux in those offset spec- 
tra is roughly half that in the center, consistent 
with quite compact emission in the central region. 

A small map of the CO J=3-2 emission is shown 
in Figure 4. For the central position we derive an 
integrated intensity of 84 K km s _1 , which for a 
conversion factor of 35.5 Jy K _1 (T^) corresponds 
to a flux of 2980 Jy km s _1 . When we convolve 
this map to simulate the 21" beam of the CO 
J=2-l spectrum, we derive an integrated intensity 
of 45 K km s" 1 or 3700 Jy km s _1 . This into- 
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Fig. 4. — Map of CO J=3-2 emission with a resolu- 
tion of 15". The integrated intensity in the central 
spectrum is 84 K km s _1 or 2980 Jy km s _1 

grated intensity is somewhat larger than the value 
of 32 K km s _1 obtained in a similar beam by 
Gerin & Phillips (1998), but agrees well with the 
measurement of Mauersberger et al. (1999). Com- 
bining the convolved spectrum with the CO J=2-l 
data gives a CO J=3-2/J=2-l ratio of 0.85 ± 0.24, 
where the uncertainty includes the estimated 20% 
calibration uncertainty in each line. 
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Fig. 5.— HCN J=4-3 with a resolution of 15". 
The integrated intensity is 7.1 K km s _1 or 
260 Jy km s" 1 . 

The HCN spectrum is shown in Figure 5, 
with an integrated intensity of 7.1 K km s _1 
(260 Jy km s" 1 ) for the HCN J=4-3 transition. 
We can combine the HCN and CO data obtained 
with similar beams to obtain an HCN J=4-3/CO 
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J=3-2 line ratio of 0.085 ± 0.024. This line ra- 
tio is quite similar to the value of 0.12 obtained 
by Solomon, Downes, & Radford (1992) for the 
J=l-0 transitions in 23-28" beams. 

3. Analysis of Single Baseline Interfero- 
metric Data 

The CSO-JCMT interferometer was used to 
measure the cross-correlated amplitude and phase 
of Arp 220 as a function of observing frequency 
and as a function of (u, v) position (Figure 1), 
which is related to the projected baseline. With a 
fixed single baseline interferometer, we cannot ob- 
tain sufficient coverage of the (u, v) plane to pro- 
duce an image of the source on the sky. Never- 
theless, the (u, v) data themselves contain a lot of 
information, which can be extracted by compar- 
ing simple models to the data (such as a single ex- 
tended disk, two point sources, etc.). We analyzed 
the amplitude and phase data separately following 
the guidelines in Lay (1994b). The method of fit- 
ting models to the visibility phases is explained in 
§ 3.1, that of fitting models to the visibility am- 
plitudes in § 3.2 and the best model fits to our 
continuum and CO J=3-2 spectral line data are 
presented in § 3.3. Readers who are interested 
primarily in the results of our analysis may wish 
to skip directly to § 3.3. 

3.1. Fitting Models to Visibility Phase 

The interferometer measures the difference in 
path lengths from the source to the two antennas 
in terms of phase. When corrected for geometric 
and instrumental effects, this phase contains infor- 
mation about the source's position in the sky. Un- 
fortunately, the CSO-JCMT interferometer is not 
stable enough to determine accurately an absolute 
position of a source from its phase. However, with 
our spectral line data, it is possible to use the red- 
shifted emission of Arp 220 as a phase reference 
and to measure the phase of the blueshifted emis- 
sion relative to the redshifted emission. Fitting 
models to this phase difference as a function of 
(u, v) position can determine the position offset 
on the sky of the blueshifted emission relative to 
the redshifted emission. As an example, for a pure 
east- west baseline, a position offset in right ascen- 
sion (RA) will result in a phase change propor- 
tional to the cosine of the hour angle of Arp 220, 



whereas an offset in declination introduces a phase 
change proportional to sin(RA) sin((5), where 5 
is the declination of Arp 220. More accurately, 
A(f> = 2w(— uARA — vAS), where u and v are the 
(u, v) coordinates of the phase reference source, 
which depends on all three components of the 
baseline as well as the declination and hour angle 
of the phase reference source (Thompson, Moran, 
& Swcnson 2001). 

3.2. Fitting Models to Visibility Ampli- 
tudes 

The visibility amplitude changes as a function 
of time depending on the structure of the source 
and the length and direction of the projected base- 
line. The visibility of a point source, i.e. a quasar, 
will be constant after correction for instrumen- 
tal and atmospheric effects. Therefore the visi- 
bilities of a point source, here 3C 345, are used 
to calibrate out any instrumental or atmospheric 
changes. An extended source might be partially 
resolved out around transit, when the projected 
baseline is the largest, in which case its visibil- 
ity amplitude would decrease. The visibility am- 
plitude of two point sources will be a minimum 
whenever their separation is a multiple of half the 
fringe spacing. From these very simple consider- 
ation, it is already clear from Figure 6 that our 
data are neither consistent with the visibility ex- 
pected from a single point source nor are the data 
in panels (a), (d) and (e) consistent with that from 
a single extended source. 

Therefore we attempted to fit the Arp 220 data 
primarily with models consisting of two sources. 
Each source is assumed to have a Gaussian bright- 
ness distribution of an elliptical shape; this ap- 
proach results in a model with 10 parameters (the 
separation of the two sources and the position an- 
gle of this separation; the flux, minor axis, major 
axis, and position angle of each source). It is pos- 
sible that models with even more sources would 
fit the data even better. However, with measure- 
ments of the visibility amplitude at only 6 or 7 
different (u, v) positions, we do not have enough 
data to constrain models with many more free pa- 
rameters. 

We used the Maximum Likelihood method to 
determine which source model fit the data best. 
We first calculated the expected visibility ampli- 
tude for each model. Since the noise of the data 
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Fig. 6.— Visibilities and model fits to the Arp 220 
data. The data are represented by squares, the 
model is the thick line and the dots are simulated 
data given the model and the measured noise, (a) 
Continuum data at 352.5 GHz. (For clarity the 
y-axis of the continuum data only runs from to 
0.6 Jy in contrast to to 2 Jy in all other panels of 
the line data.) The model represents the visibility 
of a binary with two point-like sources with fluxes 
of 0.25 Jy and 0.15 Jy, and separated by l'.'O at a 
position angle east of north (P.A.) of 80°. (b) CO 
J=3-2 visibilities (after continuum subtraction) in- 
tegrated over 5120 to 5280 km s _1 (radio defini- 
tion) and the binary model fit to the data. (Source 
1: 3.2 Jy, 2'.'4 x 274; Source 2: 0.83 Jy, 1" x 1"; 
Separation: 2'.'8 at P.A. 30°). (c) CO J=3-2 vis- 
ibilities integrated over 5280 to 5440 km s _1 and 
binary model fit to the data. (Source 1: 0.65 Jy, 
079x079; Source 2: 1.0 Jy, 0795 x 0795; Separation: 
170 at P.A. 52°). (d) CO J=3-2 visibilities inte- 
grated over 5440 to 5600 km s _1 and binary model 
fit to the data. (Source 1: 0.55 Jy, 0784 x 0784; 
Source 2: 0.92 Jy, 170 x l'.'O; Separation: 1745 at 
P.A. 115°). (c) CO J=3-2 visibilities integrated 
over 5600 to 5720 km s" 1 and binary model fit to 
the data. (Source 1: 0.3 Jy, point-like; Source 2: 
0.1 Jy, 076 x 076; Separation: 177 at P.A. 130°). 



is known, the probability of each data point given 
the model can be calculated. In the case of high 
signal to noise, the probability is described by 
a Gaussian distribution. However, in our case, 
the noise is comparable to the signal and there- 
fore the Rice distribution describes our data better 
(Thompson et al. 2001). The product of the prob- 
ability of each of the data points given the model 
gives the likelihood of the model. The model with 
the greatest product of probabilities is the most 
likely. 

For the continuum data we then calculated the 
Bayesian error of each of the ten parameters for 
the most likely model. First, the probabilities for 
each model were normalized by dividing them by 
the maximum probability. Then the probabilities 
of all models with the parameter a having a given 
value oi were summed. (This sum is the marginal 
probability of a\, which is independent of the other 
nine parameters.) This process was repeated for 
all the different values 02, a 3 etc. of a. (If we had 
a continuous rather than a discrete set of values 
for a, this process would be equivalent to integrat- 
ing the probability p(a, b, c, ...,j) of the ten dimen- 
sional parameter space over all parameters but a 
to give the one dimensional marginal probability 
of a, p mar (a) = J...JJp(a,b,c,...,j)dbdc...dj.) 
Usually, but not necessarily, the value a max with 
the highest marginal probability is also the value 
found for the most likely model. 

In a small region around a maxi the probabil- 
ity distribution can be approximated by a Gaus- 
sian. We defined the one sigma Bayesian error 
of o-max as the distance between a max and a a , 
where a a is the location where the marginal prob- 
ability drops to ex_p(-l/2)p mar (a max )- (More 
precisely, a„ should define the boundaries within 
which 68% of the marginal probability lies. In the 
case where the marginal probability can be ap- 
proximated by a Gaussian, p m ar{a<r) will have the 
value exp(— l/2)p mar (a max ). Since in our case the 
marginal probability distribution of the continuum 
data appears Gaussian, we used the latter method 
to determine the error.) The error on each param- 
eter gives an indication of how well the measured 
data determine each parameter in a given model, 
but the error contains no information about the 
validity of the model. 

For the line data the Bayesian error could not 
be calculated, because the marginal probability is 
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not Gaussian, as there are large secondary max- 
ima of the probability. We have therefore taken 
a different approach: after finding the most likely 
model, we fixed all parameters but one and defined 
the la error as the distance where the probabil- 
ity falls to cxp(— l/2)p(a max ). The difference to 
the continuum method is that we do not integrate 
over all other parameters, i.e. we are assuming 
that all other parameters are accurately known. 
The advantage of this method is that it ignores 
most secondary maxima, as we are only taking a 
one-dimensional cut in the 11-dimensional prob- 
ability surface. Since we are not taking account 
of all other planes, the errors are much smaller 
and do not reflect the error due to the existence 
of secondary maxima. (It might help to picture 
a two-dimensional probability distribution as a 
landscape with mountains, where the two dimen- 
sional plane represents the two parameters and the 
height the probability. The case of the continuum 
data can be imagined as one very high mountain 
in a landscape of low hills; the line data would 
correspond to a complicated mountain range with 
many nearly equally high mountains.) 

From the maximum likelihood method, we can 
only determine which of the models we tried is the 
most likely. It is, however, possible that an en- 
tirely different model we have not considered fits 
the data much better. For this reason it is very 
helpful to have maps of the source at other wave- 
lengths (e.g. around 230 GHz), which show the 
morphology of the source. 

3.3. Results of Model Fitting to Contin- 
uum and CO 3—3-2 Data 

The 342.5 GHz continuum measurements of 
Arp 220 were fit using both single and binary 
source models. A single source model did not 
fit the data well. Figure 6a shows the visibility 
data at 342.5 GHz and the best binary fit to the 
data. The parameters of the best binary fit are 
listed in Table 3.3, which also gives a comparison 
to published results from continuum data at other 
wavelengths. Note in particular that, while the 
continuum data can constrain the strength and 
separation of the two sources, the models cannot 
tell us whether the eastern or the western source 
is the stronger one. This limitation is due to the 
lack of absolute phase information as discussed in 
§ 3.1. 
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Fig. 7. — Measurement of source separa- 
tion using the phase difference between the 
blueshifted (5290-5440 km s" 1 ) and redshifted 
(5550-5810 km s" 1 ) emission of Arp 220. The 
top graph shows the phase of the redshifted emis- 
sion, the bottom graph shows the phase of the 
blueshifted emission. In both cases, a line fitting 
the phase of the redshifted emission has been sub- 
tracted; the scatter in the red phase therefore rep- 
resents the noise in the data. The model fit is rep- 
resented by the thick line and suggests an offset 
of (0."08,-0."12) for the redshifted emission from 
an arbitrary phase reference and (— 1."28, 0."18) 
for the blueshifted emission. 
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Table 1 
Continuum Data 





These data 


S99 


D&S 


Scoville 


Baan & Haschick 




342.5 GHz 


230 GHz 


230 GHz 


2.2^im 


4.83 GHz 


Morphology 


binary 


binary 


binary 


three sources 


binary 


Total Flux (mjy) 


400±40 


208 


175 




214 


Separation (arcscc) 


1.0 ± 0.1 


0.9 


0.8 


1.13 (NE-W) 
1.05 (SE-W) 


0.98 


Position 


80 ± 12 


100 


100 


86 (NE-W) 


98 


Angle (deg) 








108 (SE-W) 




Eastern Source: 












Flux mjy 


150±10 


66 


30 




88 


Size (arcsec 2 ) 


< 0.6 2 


< 0.2 2 


0.6 x 0.6 


0.34 x 0.26 (NE) 
0.26 x 0.23 (SE) 


0.42 x 0.31 


Western Source: 












Flux (mjy) 


250±40 


142 


90 




112 


Size (arcsec 2 ) 


< 0.6 2 


0.32 x 0.19 


0.3 x 0.3 


0.49 x 0.22 


0.30 x 0.21 


Disk 












Flux (mjy) 
Size (arcsec 2 ) 








55 




14 












Flux cal. error 


20% 


20% 


20% 







Note. — Parameters of the binary model that fits the 342.5 GHz continuum data best. The errors quoted 
are the 3cr Baycsian error of the model fitting (see § 3.2). The fluxes have an additional 20% calibration 
error. For comparison measurements at millimeter (S99; D&S), near-infrared (Scoville et al. 1998), and 
radio wavelengths (Baan & Haschick 1995) are also listed. The radio fluxes quoted are from the naturally 
weighted 4.83 GHz continuum maps. 



To determine the relative location of the two 
CO emission peaks, we divided the CO J=3-2 
spectra into two velocity bins. We chose a velocity 
range (radio definition) from 5290 - 5440 km s _1 
for the blue component, and 5550 - 5810 km s _1 
for the red component. We fit a smooth "gain" 
curve to the phase of the red emission versus time 
and then subtracted this curve from both the red 
and blue phase. The results of this phase refer- 
encing are shown in Figure 7. We then fit mod- 
els to both the red and the blue phase data us- 
ing the method described in § 3.1, with the result 
that the averaged blue emission originates from 
a location offset by — 1."36 ± 0."14 in right as- 
cension and 0."30 ± 0."14 in declination from the 
averaged red emission. In the simple case where 
all of the redshifted emission is emitted from one 
point-like component and all the blueshifted emis- 
sion from another, this result gives the separation 
of the two components. However, the analysis 
of the CO J=3-2 amplitude data described below 
shows that this model of two point sources may be 
an over-simplification. Therefore, the result from 
the phase data may only give a rough estimate 
of the separation of the two locations producing 



most of the CO J=3-2 emission. From the phase 
data, it is certain, however, that the redshifted and 
blueshifted CO J=3-2 emission originate from dif- 
ferent locations in Arp 220. 

Published lower frequency maps show CO J=2- 
1 and J=l-0 emission with a complex morphology 
(D&S; S99), which may be difficult to model with 
simple single or binary source models. To simplify 
the morphology to some extent, we divided the CO 
data in four velocity bins of 5120 - 5280 km s _1 , 
5280 - 5440 km s~\ 5440 - 5600 km s" 1 and 5600 - 
5760 km s , each of which corresponds to four ve- 
locity channels in Figure 20 from D&S. But even 
in these narrower velocities bins, the CO J=2-l 
emission is complex. In addition, our single (u, v) 
track cannot sample the more extended emission 
adequately, nor does it have much north-south res- 
olution. 

The visibility amplitude plots themselves show 
some clear signatures of the source structure. Fig- 
ure 6b-d all show an increase in amplitude around 
an hour angle of 5 h, which indicates extended 
emission in the east-west direction that was re- 
solved out at longer baselines. In contrast, the 
highly redshifted component of the line in Fig- 



9 



lire 6e does not show a rise in amplitude as Arp 220 
sets, which means that there is no extended (arc- 
second scale) east-west component at these veloci- 
ties. The continuum data in Figure 6a also do not 
show a large increase in amplitude at this hour 
angle. 

For the line data (in contrast to the contin- 
uum data), there are many very different models 
which fit each set of visibility amplitudes fairly 
well. Of these models, we have chosen the one 
which seemed most consistent with the CO J=2-l 
maps of D&S. The visibility data and the model 
fits are shown in Figure 6 and the parameters of 
the model are listed in Table 3.3. The errors in 
Table 3.3 are 3<r errors of each parameter assum- 
ing all other parameters are accurately known (see 
§ 3.2). In addition to those errors there are: (1) 
an error due to the existence of secondary maxima 
(i.e. a completely different combination of param- 
eters might give a nearly equally good fit); (2) 
an error due to only trying binary models; (3) an 
error due to the uncertainty in attributing the vis- 
ibilities to the different features; and (4) the fluxes 
have an additional 25% calibration error. 

Each model is only one possible way of describ- 
ing our data and parameters in Table 3.3 should 
be used with caution. The visibility amplitudes 
are best fit by binaries. However, these binaries 
are not the east- west binary seen in the continuum 
data; rather, it appears to be one or the other of 
the two continuum nuclei plus a much more ex- 
tended CO disk. This larger CO disk appears at 
different separations and position angles from the 
more compact nuclei in different velocity channels, 
which complicates the interpretation of the results. 

When comparing the CO J=3-2 emission to 
maps, it is best to use the visibilities directly 
rather than the results from the model fitting. 
Thus, we extrapolated the visibilities from the CO 
J=2-l OVRO data (S99) to the (u, v) positions 
observed in CO J=3-2 with the CSO-JCMT for a 
direct comparison, which is shown in Figure 8. 

4. Discussion 

4.1. Continuum Data 

The binary model parameters that fit the con- 
tinuum data best are listed in Table 3.3, which 
also gives a comparison to continuum data at 
other frequencies from the literature. The binary 
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Fig. 8.— Comparison of CO J=3-2 and CO J=2-l 
visibilities. The visibilities were calculated for the 
(u, v) positions shown in Figure 1 and for four ve- 
locity bins of 5120-5280, 5280-5440, 5440-5600 and 
5600-5760 km s" 1 (radio definition) for CO J=3-2 
(triangles, dashed line) and CO J=2-l (hexagons, 
solid line). The calibration error on the CO J=3-2 
data is estimated to be 25%, that of the CO J=2-l 
data 10%, respectively. 
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Table 2 

Best Fitting Models to CO 3-2 Interferometer Data 



CO .7=3-2 


5120-5280 


5280-5440 


5440-5600 


5600-5760 


Vcl. Bin 


km s — 1 


km s — 1 


km s~ 1 


km s~ 1 


Separation (arcscc) 


2.8±0.3 


1.0±0.1 


1.45±0.1 


1.7±0.3 


Position Angle (deg) 


30±6 


52±9 


115±6 


130±12 


Eastern Source 










Flux (Jy) 






0.55±0.15 


0.3±0.06 


Diameter (arcscc) 






0.84±0.12 


<0.6 


Western Source 










Flux (Jy) 


0.83±0.18 


0.65±0.3 


0.92±0.12 




Diameter (arcscc) 


1.0±0.3 


0.9±0.25 


1.0±0.12 




Disk 










Flux (Jy) 


3.2±0.6 


1.0±0.15 




0.1±0.12 


Size (arcsec 2 ) 


2.4±1.2 


0.95±0.12 




0.6±0.6 


Pos. Angle (deg) 








90 



Note. — Parameters of models fitting the visibilities in four different velocity 
ranges of CO J— 3-2 emission. The data do not constrain the models well and we 
have chosen the model that is most consistent with the CO J— 2-1 maps of D&S. 
The errors quoted here are the 3(T errors of each parameter under the assumption 
that all other parameters are correct (see § 3.2). There are additional, mainly 
systematic, errors discussed in § 3.2 such that the overall error might be as large 
as a factor of 2. 



separation determined from the continuum data 
is slightly larger and the position angle slightly 
smaller compared to the two sources seen in the 
1.3 mm continuum map of S99 and D&S. Our 
data, however, agree reasonably well with the po- 
sition of the formaldehyde emission, whose peaks 
are separated by l'.'l at 79° (Baan & Haschick 
1995). Formaldehyde is associated with star burst- 
ing gas (Baan & Haschick 1995), which will also 
contain hot dust emitting at submillimeter wave- 
lengths; it is therefore not surprising that the sub- 
millimeter emission shows the same morphology 
as the formaldehyde emission. The two contin- 
uum sources are point-like and their size cannot 
be resolved with the CSO-JCMT interferometer. 
From the model fitting (§3.2), the 3a upper limit 
of 0?6 is consistent with the source sizes seen by 
D&S and S99 as well as the 2.2 ^m emission, which 
traces dust-enshrouded young stars Scoville et al. 
(1998). 

We detect only about half of the single dish 
continuum flux [extrapolated with (3—1 from 
the 1.1 ± 0.4 Jy total flux, which was measured 
by Ealcs, Wynn- Williams, & Duncan (1989) at 
375 GHz with the JCMT]. This result is slightly 
surprising because interferometry maps do not 
show significant extended continuum emission at 



230 GHz. Most of the measured 400 mJy ± 
80 mJy continuum flux is expected to come from 
dust emission rather than synchrotron or free-free 
emission. The flux ratio of the two sources is about 
1.7, similar to the 230 GHz data from S99, which 
suggests that the continuum flux arises under sim- 
ilar physical conditions in both sources. (D&S di- 
vide the flux into 3 sources and so a direct com- 
parison is more difficult.) To first order, the con- 
tinuum flux is proportional to the product of the 
dust mass and the dust temperature. Unless the 
dust temperatures vary by more than a factor of 3 
between the two nuclei, the mass of dust of the two 
nuclei must be of the same order of magnitude. 

Gas and dust masses, column densities, volume 
densities and visual extinctions are listed in Ta- 
ble 3. To calculate the dust masses we adopted a 
temperature of 42 K, which was derived by Scov- 
ille et al. (1991) by fitting a black body curve to 
IRAS data for Arp 220. We used an opacity coef- 
ficient of 1 cm 2 g _1 . The gas mass was calculated 
using a gas to dust mass ratio of 100 as suggested 
by SYB. Our continuum measurements only give 
an upper limit of the source diameter of 0'.'6. We 
used this value for both the east and west source 
to compute the lower limits of the column densi- 
ties, volume densities and visual extinctions. We 
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Table 3 

Masses, densities and extinction 



Frequency 


342 GHz 


230 GHz 


Reference 


These Data 


D&S 1 


Eastern Source 






M duet (M ) 


2.8 x 10 7 




M gas (M ) 


2.8 x 10 9 


0.6 x 10 9 


Njf 2 (cm 2 ) 


> 5.5 x 10 24 




PH 2 ( cm_3 ) 


> 15000 


900-20000 a 


Av (mag) 


>5800 




Western Source 






M dust (M s ) 


4.7 x 10 7 




M gas (M ) 
Nh 2 (cm- 2 ) 


4.7 x 10 9 


1.1 x 10 9 


> 9 x 10 24 




PH 2 (cm -3 ) 


>25000 


900-22000* 


A v (mag) 


>9600 




Disk 






Nh 2 (cm- 2 ) 




10 24 


A v (mag) 




-1000 



a Dcrivcd by fitting models of a disk of chang- 
ing thickness to CO data. 



References. — (1) Downes & Solomon 1998 



used N(H)/E(B-V)=5.8 x 10 21 from Bohlin, Sav- 
age, & Drake (1978) and adopt A„/E(B-V)=3.1 
to derive the visual extinction. 

The total dust mass derived for our two sources 
of 7.5 x 10 7 M Q agrees fairly well with the esti- 
mate of 5 x 10 7 M from the 110 GHz continuum 
emission by Scoville et al. (1991). However, our 
gas mass estimates for the two sources are signif- 
icantly larger than those of D&S using 230 GHz 
data; moreover they are larger than the dynami- 
cal masses of ~ 1O 9 M for these two regions (S99). 
Since the gas mass cannot be larger than the dy- 
namical mass, one of our assumptions must be 
wrong. The assumed gas to dust ratio of 100 is 
already at the low end of ranges typically adopted 
and so reducing that ratio seems unlikely. Either 
a larger grain emissivity or a larger temperature 
would act to reduce the dust masses and hence the 
gas masses derived from them. (Using 100 K ver- 
sus 42 K will reduce our mass estimates by a factor 
of 2.6.) Changes in grain properties and elevated 
temperatures would not be unexpected given the 
unusual properties of this starburst region. 



4.2. CO Data 

The single dish CO J=2-l and CO J=3-2 data 
are shown in Figures 3 and 4. The lines shapes are 
quite similar for these two transitions and the mea- 
sured line ratio of 0.85 agrees quite well with ob- 
servations of this transition in normal spiral galax- 
ies and probably indicates moderately warm (30- 
50 K) gas (Wilson, Walker, & Thornley 1997). 
The CO J=3-2 emission of Arp 220 is clearly mod- 
erately extended, as the flux within a 22" beam is 
20% larger than the flux within a 15" beam. 

Our CO J=3-2 integrated intensity agrees quite 
well with that measured by Mauersberger et al. 
(1999) in a similar beam, once both measurements 
are converted to the same temperature scale (i.e. 
both in Ta* or Tmb)- There are slight differ- 
ences in the line shape between the two observa- 
tions, particularly in the relative strength of the 
red and blue peaks; these differences are likely due 
to slightly different pointing between the two ob- 
servations, as the rms pointing accuracy in the 
Mauersberger et al. data is 5". 

We can combine our CO data with the CO J=l- 
measurement from Solomon, Radford, & Downes 
(1990) to obtain beam-matched estimates of the 
CO J=3-2/2-l and J=2-l/l-0 line ratios. These 
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Fig. 9. — Large velocity gradient models for Arp 
220 for four different kinetic temperatures. Three 
lines are shown for the 12 CO J=3-2/2-l and J=2- 
1/1-0 line ratios, which correspond to values of 
0.84 ± 0.24 and 0.82 ± 0.23, respectively. For 
the 12 CO/ 13 CO J=l-0 line ratio, two lines are 
shown, which correspond to the measured 3c up- 
per limit of 19 from Aalto et al. (1991) and this 
value minus 30% to account for calibration un- 
certainty. For this line ratio, an arrow indicates 
the direction of the allowed region. For a given 
temperature, allowed values of density and col- 
umn density occur where the three allowed re- 
gions for the three line ratios intersect. For ex- 
ample, for Tk = 100 K, one allowed region is 
the triangular region with nn 2 ~ 10 3 cm" 3 and 
N( 12 CO)/dv ~ 10 17 cm" 2 km" 1 s, and there is a 
second allowed region with n H2 > 10 4 cm" 3 and 
N( 12 CO)/dv ~ 10 18 cm" 2 km" 1 s. 



two line ratios alone do not allow us to place useful 
constraints on the physical conditions in the gas 
using large velocity gradient models. However, we 
can get some interesting constraints, if we include 
the upper limit to the 12 CO/ 13 CO J=l-0 line ratio 
of 19 from Aalto et al. (1991). One caveat is that 
this isotopic line ratio was measured in a much 
larger beam (54") and so may trace emission from 
beyond the nuclear region. We compared these 
three line ratios to the output from large velocity 
gradient models (Figure 9) run with a range of 
temperatures, densities, and CO column densities 
(10-300 K, 10-10 6 cm" 3 , 10 15 -10 20 cm" 2 km" 1 s). 
Solutions could be found for all the temperatures 
we investigated; in general, we can place only a 
lower limit on the density for a given temperature, 
and the minimum allowable density decreases as 
the temperature increases. 

The interferometric spectra are shown in Fig- 
ure 2 and the visibilities in Figure 6. Figure 2 
shows that the shapes of the spectra change dras- 
tically with (u, v) position. In particular, strong, 
mainly blucshiftcd emission can be seen in the last 
spectrum; since the fringe spacing is 5. "8, this 
emission must be quite extended. In addition, the 
visibility phases clearly indicate that there are at 
least two sources of emission (Figure 7). 

The models fitting the visibilities best are listed 
in Tables 3.3 and 4. However, the values crucially 
depend on the model chosen and should be used 
with caution. The tabulated models suggest a to- 
tal flux of 1210 Jy km s" 1 , which is 40% of the 
3000 Jy km s" 1 observed with the JCMT in a 15" 
beam. A lower limit on the flux can be obtained 
by simply adding the visibilities in different ve- 
locity bins observed for the 5"8 spacing which re- 
sults in 700 Jy km s _1 . The comparison of single 
dish and interferometer data demonstrates that a 
significant amount of the CO J=3-2 flux is more 
extended than can be detected with a 5"8 maxi- 
mum fringe spacing, and care should be taken to 
compare only data with the same beam size/fringe 
spacing. 

Using the empirical relationship between CO 
emission and molecular mass (Petitpas & Wil- 
son 1998), the average CO J=3-2/l-0 line ratio 
of 0.7 and the fluxes determined from the model 
fitting (Table 4) we obtained molecular masses of 
5.3 x 10 9 M for the western source, 1.9 x 10 9 M 
for the eastern source and 9.6 x 10 9 M Q for the 
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Table 4 

Comparison of Interferometric CO Observations 



CO Transition 


CO 3-2 


CO 2-1 


CO 1-0 


Reference 


These Data 


D&S 1 


D&S 1 


Separation (arcscc) 


1.45±0.1 


1.3±0.1 




Position Angle (deg) 


115±6 


85±6 




Eastern Source 








Flux (Jy km s" 1 ) 


140±30 


220±44 




Source diameter (arcsec) 


0.54±0.23 


0.9±0.1 




Western Source 








Flux (Jy km s _1 ) 


380±60 


130±26 




Source diameter (arcsec) 


0.97±0.13 


0.3±0.1 




Disk 








Flux (Jy km s _1 ) 


690±90 


750±150 




Source diameter (arcsec) 


2.0±0.9 


1.8±0.1 




Total 








Flux (Jy km s -1 ) 


1210±110 


1100±220 


410±82 



Note. — Averaged parameters of models from Table 3.3. (The average 
source sizes are not the simple arithmetic average but are weighted by 
fluxes.) The errors quoted for our data are the 3(T errors from Table 3.3 
propagated according to Gaussian error analysis. The fluxes have a 25% 
calibration error. In addition, there arc systematic errors discussed in 
§ 3.2. The errors of the CO J=2-l and J=l-0 data arc taken from D&S 
and are expected to represent the overall error. 



References. — (1) Downes & Solomon 1998 



disk. These masses agree well with the masses de- 
rived from the continuum emission in § 4.1, but are 
still larger than the dynamical masses; this result 
suggests the CO-to-H2 conversion factor is lower 
than the value in the Milky Way consistent with 
other results for ultraluminous infrared galaxies 
(Solomon et al. 1997). 

The CO J=3-2 interferometric data suggest 
that the western (blucshifted) source is brighter 
than the eastern (redshifted) source; the sin- 
gle dish spectrum also indicates slightly stronger 
emission in the blucshifted part of the line. This 
interferometric result is in the opposite sense to 
that seen by D&S, and could indicated tempera- 
ture differences between the two nuclei. [However, 
S99 clearly show the western source is brighter 
than the eastern source at CO J=2-l, so perhaps 
there is a typographical error in D&S.] 

Morphologically, the CO J=3-2 interferometric 
data indicate the presence of two fairly compact 
emission regions with a more extended disk. How- 
ever, the morphology appears to be complex. In 
general, our data are consistent with D&S as well 
as S99 data and small inconsistencies are most 



likely due to an over-simplification in our model, 
which fits only two Gaussian peaks to the emis- 
sion. 

Figure 8 shows visibilities of CO J=3-2 and 
J=2-l taken at exactly the same (u, v) positions. 
The shape of the visibility curves roughly agree 
in the two most redshifted bins, with larger dif- 
ferences between the two transitions seen in the 
two blueshifted bins. In particular, in the 5280- 
5440 km s _1 frequency bin, we see more extended 
emission in CO J=2-l than CO J=3-2. The CO 
J =2-1 might also show more evidence for a binary 
structure. The visibilities at HA=-0.8 are the best 
representatives of the eastern and western source 
we have, as the disk will be mostly resolved out at 
this long projected baseline (-170kA,60kA, fringe 
spacing l'.'l). For this (u, v) point the CO J=3- 
2/J=2-l line ratios (converted to K scale) are ap- 
proximately 0.14, 0.12, 0.08 and 0.92 for increasing 
velocities. In the high resolution CO J=2-l map of 
S99, the western source dominates between 5050 
and 5450 km s _1 , which corresponds to our first 
two velocity bins, and the eastern source domi- 
nates between 5500 and 5650 km s^ 1 , which cor- 
responds to our last two velocity bins. The average 
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line ratio for the blucshiftcd western source is 0.1, 
significantly below the average of 0.5 for the east- 
ern source. This results could be interpreted as the 
eastern source being warmer (or denser) than the 
western source and is consistent with results from 
single dish HCN and CO observations discussed in 
the next section. 

4.3. HCN Data 

Figure 5 shows the HCN J=4-3 spectrum ob- 
tained with the JCMT. In contrast to the CO 
J=3-2 spectrum (Figure 3) and the HCN J=l-0 
spectrum (Solomon et al. 1992), the HCN J=4-3 
spectrum is dominated by a single rcdshifted emis- 
sion peak. We divided the HCN emission into red 
and blueshifted emission at 5430 km s _1 , where 
the CO as well as HCN 1-0 spectra have a dip 
in the intensity. We converted the HCN J=l-0 
line to the temperature scale using a main 
beam efficiency of 0.6 (Radford et al. 1991) and 
scaled the integrated intensity by a factor of four 
to correct for the difference in beam sizes. Note 
that the difference in beam size leads to an uncer- 
tainty in the line ratio because we do not know the 
true source structure; the factor of four scaling is 
only strictly appropriate for a point source. Our 
HCN J=4-3 and CO J=3-2 data were obtained 
with almost identical beam sizes and so we can 
calculate a line ratio without the need to apply 
any additional corrections. The integrated inten- 
sities and the line ratios of HCN J=4-3/J=l-0 and 
HCN J=4-3/CO J=3-2 are listed in Table 5. 

The line ratios both in HCN J=4-3/J=l-0 and 
HCN J=4-3/CO J=3-2 are larger for the red- 
shifted eastern source than the blueshifted west- 
ern source. Since the HCN J=4-3/J=l-0 line ra- 
tio is not affected by abundance changes, it seems 
likely that the physical conditions are different be- 
tween the two emission peaks and that the eastern 
source is either at a higher temperature and/or has 
a higher density. These conclusions are supported 
by the work of (Aalto et al. 2002), who detected 
CN J =2-1 emission from the blueshifted part of 
the line, and HC3N emission in the rcdshifted part 
of the line. They suggest that CN emission is an 
indicator of a photon dominated region while the 
HC3N emission is an indicator of hot cores, and 
suggest that the two nuclei may be in different evo- 
lutionary states. Infrared maps show a dust lane 
across the eastern nucleus and a more symmetric 



morphology in the western nucleus (Scoville et al. 
1998), which could be at least superficially consis- 
tent with our observed line ratio variations. D&S 
suggest that the western nucleus is undergoing an 
intense starburst; if this starburst has dispersed 
much of the gas in the western nucleus, while the 
eastern nucleus continues to contain a larger mass 
of gas at high densities, this scenario would also 
be consistent with our data. 

5. Conclusions 

We have presented the first interferometric 
observations of Arp 220 at submillimeter wave- 
lengths. The interferometric visibilities of the CO 
J=3-2 line and 342 GHz continuum are largely 
consistent with the emission morphology seen pre- 
viously at lower frequencies. We clearly detect 
continuum and CO J=3-2 emission from at least 
two sources separated by ~ 1" at P. A. ~ 80°. The 
CO J=3-2 visibility amplitudes show additional 
extended structure with a complex morphology. 
Masses, column densities, volume densities and 
optical extinction calculated for both emission 
sources agree with previous estimates within the 
errors and underline that the center of Arp 220 
contains large amounts of molecular gas (~ 8 10 9 
M ). 

Single-dish data indicate that the CO J=3-2 
emission is moderately extended compared to the 
15" beam of the JCMT. Though the continuum 
visibilities show no sign of an extended source, 
the single dish continuum flux is about twice that 
detected with the interferometer. In single dish 
data, the HCN J=4-3/J=l-0, HCN J=4-3/CO 
J=3-2 and CO J=3-2/J=2-l ratios are all larger 
for the redshifted portion of the line than for the 
blueshifted portion. These observations suggest 
that the redshifted eastern source is denser and/or 
hotter than the blueshifted western source. This 
results could provide support for D&S hypothesis 
that the western source is currently undergoing 
an intense starburst that has dispersed the dense 
gas, whereas the eastern source still harbors dense 
molecular material. 

We are extremely grateful to Kazushi Sakamoto 
who extracted a CO J=2-l visibility from his data 
that exactly matched our (u, v) track. Further- 
more we would like to thank David Brown, Claire 
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Table 5 

hcn intensities and line ratios 



Source 


HCN J=4-3 
(K km s" 1 ) 


HCN J=l-0 a 
(K km s- 1 ) 


CO J=3-2 
(K km s" 1 ) 


HCN J=4-3/J=l-0 


HCN J=4-3/CO J=3-2 


Eastern Source 


4.0 


9.6 


36 


0.42 


0.11 


Western Source 


3.0 


10.8 


48 


0.28 


0.064 


Total 


7.1 


20.4 


84 


0.35 


0.085 



a Thc integrated Tmb from Solomon ct al. (1992) were converted into Ta* by multiplication with the telescope 
efficiency of 0.6 and corrected to a 14" beam by multiplying by a factor of four. 
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ported by NSF grant AST 99-80846. The JCMT 
is operated by the Royal Observatories on behalf 
of the Particle Physics and Astronomy Research 
Council of the United Kingdom, the Netherlands 
Organization of Scientific Research, and the Na- 
tional Research Council of Canada. MCW was 
supported by an HSP III scholar-ship of the Ger- 
man Academic Exchange Service as well as the Sir 
Isaac Newton Scholarship. 

REFERENCES 

Aalto, S., Johansson, L. E. B., Booth, R. S., & 
Black, J. H., 1991, A&A, 249, 323 

Aalto, S., Polatidis, A. C, Hiittemeister, S., & 
Curran, S. J., 2002, A&A, 381, 783 

Arp, H. 1966, Atlas of Peculiar Galaxies 
(Pasadena: Caltcch) 

Baan, W. A., & Haschick, A. D. 1995, ApJ, 454, 
745 

Bccklin, E. E., & Wynn- Williams, C. G. 1987, in 
Star Formation in Galaxies, ed. C. J. Lonsdale 
(Washington: NASA), 643 

Bohlin, R. C, Savage, B. D., & Drake, J. F., 1978, 
ApJ, 224, 132 

Condon, J. J., Huang, Z. P., Yin, Q. F., & Thuan, 
T. X. 1991, ApJ, 378, 65 

Downes, D. & Solomon, P. M. 1998, ApJ, 507, 615 

Eales, S. A., Wynn- Williams, C. G., & Duncan, 
W. D. 1989, ApJ, 339, 859 

Eckart, A., & Downes, D. 2001, ApJ, 551,730 



Emerson, J. P., et al. 1984, Nature, 311, 237 

Genzel, R., et al., 1998, ApJ, 498, 579 

Gerin, M. & Phillips, T. G. 1998, ApJ, 509, L17 

Heckman T., Dahlcm M., Eales S. & Weaver K. 
1996, ApJ, 457, 616 

Ivison, R. J., Smail, I., Bargcr, A. J., Kncib, J.- 
P., Blain, A. W., Owen, F. N., Kerr, T. H., & 
Cowie, L. L. 2000, MNRAS, 315, 2091 

Joseph R. & Wright G. 1995, MNRAS, 214, 87 

Lay, O. P., Carlstrom, J. E., Hills, R. E., & 
Phillips, T. 1994a, ApJ, 434, 75 

Lay, O.P. 1994b, PhD thesis, Cambridge Univer- 
sity 

Lay, O. P., Carlstrom, J. E., & Hills, R. E. 1995, 
ApJ, 452, 73 

Lay, O. P., Carlstrom, J. E., & Hills, R. E. 1997, 
ApJ, 489, 917 

Lutz, D., ct al. 1996, A&A, 315, L137 

Mauersberger, R., Hcnkcl, C, Walsh, W., & 
Schulz, A., 1999, A&A, 341, 256 

Mirabel, I. F., & Sanders, D. B. 1989, ApJ, 340, 
L53 

Norris, R. P. 1988, MNRAS, 230, 345 

Padman R. 1993, 'SPECX Users Manual', 
Cavendish Laboratory, Cambridge 

Petitpas, G. R. & Wilson, C. D. 1998, ApJ, 496, 
226 



16 



Rickc, G. H., Cutri, R. M., Black, J. H., Kailey, W. 
F., McAlary, C. W., Lebofsky, M. J., & Elston, 
R. 1985, ApJ,290, 116 

Radford S., Solomon P. & Downes D. 1991, ApJ, 
368, L15 

Robson, E. I., Stevens, J. A., & Jenncss, T. 2001, 
MNRAS, 327, 751 

Sanders, D. B., Scovillc, N. Z., & Soifcr, B. T. 

1991, ApJ, 370, 158 

Sakamoto, K., Scoville, N. Z., Yun, M. S., Crosas, 
M., Genzcl, R., & Tacconi, L. J. 1999, ApJ, 
514, 68 

Scoville, N. Z., Sargent, A. I., Sanders, D. B. , & 
Soifcr, B. T. 1991, ApJ, 366, L5 

Scoville, N. Z., Yun, M. S., & Bryant, P. M. 1997, 
ApJ, 484, 702 

Scoville, N., Evans, A., Thompson, R., Riekc, M., 
Schneider, G., Low, F., Hines, D. & Stobbie, B. 
1998, ApJ, 492, L107 

Smith H., Lonsdale C, Lonsdale C. & Diamond 
P. 1998, ApJ, 493, L17 

Soifcr, B. T., et al., 1987, ApJ, 320, 238 

Sopp, H. M. & Alexander, P. 1991, MNRAS, 251, 
112 

Solomon, P. M., Radford, S. J. E., & Downes, D., 
1990, ApJ, 348, L53 

Solomon, P. M., Downes, D., & Radford, S. J. E. 

1992, ApJ, 387, L55 

Solomon, P. M., Downes, D. Radford, S. J. E., & 
Barrett, J. W., ApJ, 478, 144 

Sturm, E., ct al., 1996, A&A, 315, L133 

Thompson, A. R., Moran, J. M., & Swenson, G. 
W. 2001, Intcrferometry and Synthesis in Ra- 
dio Astronomy (2nd ed.; New York, NY: John 
Wiley & Sons) 

Wilson, C. D., Walker, C. E., Thornlcy, M. D. 
1997, ApJ, 483, 210 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.0. 



